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Abstract 

One of the most powerful methods for obtaining a precision measurement of 
Mw at LEP2 is from direct reconstruction of the invariant mass distribution of 
W bosons produced in e + e~ — ► W + W~ . We investigate the effects on the W 
line shape, and in particular on the average invariant mass, of phase-space and 
first-order Coulomb corrections. The latter are shown to have a non-negligible 
effect on the reconstructed mass, inducing shifts of order —20 MeV at collision 
energies above threshold, compared to the Born approximation cross section. 
The sign and magnitude of the effect can be understood in a simple model 
calculation in which one of the W bosons is assumed to be stable. 



1 Introduction 



A precision measurement of the mass of the W boson, My/, is one of the main objec- 
tives of the LEP2 physics programme. There has already been considerable progress 
in the precision determination of My/ from W — > Iv events at the Tevatron pp collider 
U, and the level of precision is expected to increase further in the coming years @. 
It is important, therefore, that the potentiality of the My/ measurements at LEP2 
should be exploited to the full. 

An obvious requirement for the success of these precise studies is a high level of 
reliability of the theoretical predictions for the various experimental observables relat- 
ing to the different methods for measuring My/ from the process e + e~ — > Vr+W 7- !] 
This in turn requires a detailed understanding of the physical phenomena which de- 
scribe the production and decay of W bosons at LEP2, in particular of the effects 
which arise from the large W boson decay width Ty/ (~ 2.1 GeV) ||. The insta- 
bility of the W bosons can, in principle, strongly modify the standard 'stable W 
results. For example, an important role can be played by the QED and QCD ra- 
diative interferences (both virtual and real) which interconnect the production and 
decay stages. Particular attention should be paid here to the virtual contributions 
corresponding to the 'charged' particle poles for which, in contrast to the photon and 
gluon poles, there is no cancellation from real emissions || [7| (see also [R W). The 
level of suppression of the width-induced radiative effects depends on the 'degree of 
inclusiveness' of the process. Thus, in the case of the total inclusive cross section for 
e + e~ —>■ W + W~ —>■ 4/, the final-state interaction effects are cancelled up to terms of 
relative order aTy/ /My/ or a 2 s Yy/ / My/- The only exception is the contribution arising 
from the Coulomb interaction between the slowly moving W bosons [TT, [TT| (see also 



12]). If the W bosons were stable particles, the effect of the Coulomb interaction on 



the total cross section has been known for a long time [13]. The modifications to the 
Coulomb corrections which arise from the instability of the W bosons are particularly 
significant near the W + W~ production threshold, yfs ~ 2M W , but become negligi- 
ble for collision energies which satisfy \fs — 2My/ ^> r^. As explained in detail in 
Ref . , the H^-boson virtuality drastically changes the on-shell value of the Coulomb 
correction even at y's — 2My/ 3> r^, but after integration over the invariant masses 
of the two W bosons the £ stable-W result is restored far above the threshold region. 
The electroweak radiative corrections to the e + e~ — > W + W~~ total cross section 



for the production of stable W bosons are now known with high accuracy [14]. The 



instability effects are well under control throughout the energy range accessible at 
LEP2 (y/s £ 200 GeV). This includes the so-called ' colour- reconnection effects' [jl5| , |l6l, 
|TT[ - non-perturbative hadronization dynamics which may affect the W + W~ — > 4q 
decay channel and which have recently attracted much attention. 

1 A summary of these methods can be found, for example, in Refs. p, m. 
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Two different methods have been proposed for a precision measurement of My/ at 
LEP2. The first uses the method of direct kinematic reconstruction of My/ from the 
decay channels 

W + W~ -> qqqq, (1) 
W + W~ -> qqlv. (2) 

The second method uses the strong dependence of the total W + W~ cross section on 
My/ near threshold to translate a measurement of aww close to = 2Mw into 
a value of Mjy. Both methods have been, and are currently being, studied in great 



detail, see for example Refs. || [|, [18[ and references therein. Although both have their 
advantages and disadvantages, the direct reconstruction method is currently believed 
to provide the better precision on My/. Since it naturally requires a higher collision 
energy to maximise the event rate, it does not conflict with the requirements of new 
particle/physics searches. 

However, the direct reconstruction method is not without problems. For example, 
to construct the two W's from the qqqq final state in ([!]) one must in principle at- 
tribute all observed hadrons to the 'correct' parent W, a procedure which is certainly 
affected by relatively unknown QCD interconnection corrections [|16j. Since a com- 
plete description of these effects is not possible at present, one has to rely on model 
predictions rather than on exact calculations; for details see Refs. [jH], 117). Fortu- 
nately, the contribution of the interconnection effects to the systematic error on My/ 
is not expected to significantly exceed the overall systematic error (not including in- 
terconnection effects), which is currently estimated at O(30 — 40) MeV. It is an open 
question whether Bose-Einstein effects might induce a further uncertainty in the mass 



determination [IB], [19| . Such problems do not of course arise for the semi-leptonic 
channel ([|), but there the event rate is smaller and an unobserved neutrino is present. 

The direct reconstruction method is based on measurements of semi-inclusive char- 
acteristics of the final state — W boson momenta, jet-jet invariant masses, opening 
angle between jets, etc. All such quantities could well be much more sensitive to insta- 
bility effects than the total cross section. For example, as could be anticipated from 
Ref. Jill (see also below), the modification of the QED Coulomb interaction could 
induce a systematic shift in the W mass of 0(a7rTy/), which is of the same order as 
the target precision. Clearly such effects should be carefully calculated and taken into 
account in the extraction of the W mass from measured distributions. Note that since 
Coulomb forces are responsible for the QED interaction between the separate hadronic 
or leptonic final states of each W, their effects could be regarded as an example of 
'QED interconnection effects'. These effects are in fact quite universal and do not 
depend on the particular W decay channels. 

In the threshold region the Coulomb contribution dominates the instability effects, 
and the Coulomb terms can be uniquely separated from the remaining electroweak 
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corrections. At higher energies, the width-induced modifications of the differential 
distributions caused by other radiative mechanisms (for example, intermediate-final 
and final-final state radiative interferences |E|, || ) may become just as important.^ 
Moreover, it is argued in Ref. that in the relativistic domain a cancellation may 
take place between the different sources of instability effects and that, as a result, the 
stable W + W~ result may be restored. In the extreme relativistic limit, (1 — 0) <C 1, 
where (3 is the W velocity, such a cancellation appears quite natural.^ However in 
the intermediate region (3 < 1, which is relevant for the LEP2 energy range, the 
situation is less clear in our opinion, and therefore needs further detailed study. Since 
\fs — M\y <C M\y at LEP2, the non-relativistic Coulomb formulae are likely to provide 
a reasonable qualitative guide to the size of the width-induced effects. 

In this study we analyse the effect of the W boson instability, as embodied in the 



first-order Coulomb formulae of Refs. Jl0|, |TT| , on the invariant mass distribution of the 



decay products. As far as we are aware, the necessity to take Coulomb-induced distor- 
tion effects in the W mass or momentum distribution into account was first pointed 
out in Ref. Ul6| . In order to expose the direct effect of this 'QED interconnection', 
we make several simplifying assumptions in our analysis: (i) the effects of initial state 
radiation are not included, although it would be straightforward to take them into 
account using the standard techniques, (ii) we assume that the W + W~ final state can 
be fully reconstructed, and (iii) we neglect possible QCD interconnection effects. Note 
that our results also apply directly to the process 77 — > W + W~ . 

The paper is organized as follows. In Section we present a simple model calcu- 
lation in which one of the W bosons is assumed to be stable. Above threshold this 
is expected to provide a reasonable qualitative understanding. Using this model, we 
study the effect on the invariant mass distribution of the decaying W of phase-space 
and first-order Coulomb effects. In Section |^ we study numerically the realistic case 
when both W bosons are off- mass-shell. Predictions for the shift in the average mass 
with and without Coulomb corrections are presented. Finally, Section |] contains our 
conclusions. 



2 Recall that in the relativistic region the Coulomb term is neither uniquely defined nor separately 
gauge invariant. 

3 It has its origins in the conservation of 'charged' currents. Another example of the vanishing of off- 
shell effects at very high energies was discussed in Ref. |2(J (see also ) , where the gluon radiation 
pattern corresponding to top quark production and decay was discussed. For the most probable 
kinematic configuration (quasi-collinear b and t), the width dependence disappears completely when 
the effects of emission at the production and decay stages are added coherently. The same behaviour 
is expected for QED radiation off fast-moving W bosons |2^] . 
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2 A model analysis 



To elucidate the physical origin of the distortion of the W decay mass distribution 
induced by the Coulomb interaction, it is instructive to consider first a simplified 
model in which one of the W bosons is assumed to be stable and the other has the 
standard W — > ff decay modes with decay width IV0 (see also |§). 

In the non-relativistic region the differential distribution of the invariant mass- 
squared Si of the unstable W can be written as 



1 + V s ") + 20 SC 



(3) 



where (Jq(s, si, s%) is the e + e — > W + W off-shell Born cross section [24]] at centre-of- 
mass energy y/s, and p{sj) is the Breit-Wigner factor^ 



p{si) = ~ 



71 



{ Si -M^ + Si r 2 w ( Si y 



(4) 



with Tw(si) = yfs~i Tw/Mw- The coefficient 5h is the 'hard' first-order radiative 
correction, (a>/2/3)5c is the first-order Coulomb contribution (see Refs. |10, ll|), and 
(3 is the velocity of the W bosons in the centre-of-mass frame. In the non-relativistic 
approximation (3 = I'pjyfs where 



V 



s- y/s[- M W )M } 



w ■ 



(5) 



In our model calculation it will be sufficient to work in the non-relativistic approxi- 
mation. It is then straightforward to show that in this case 5c is given by the same 
formulae as in the realistic case of two unstable W bosons 
by T w /2. Thus 



10, 11], with F w replaced 



6, 



c 



7T 



2 arctan 



' \ K \ 2 ~ P 2 
2p Re(fv) 



with 



K 



'-M w (E + \iT w ). 
Here E is the non-relativistic energy of the W bosons, 



E 



2M, 



w ■ 



(6) 



(7) 



(8) 



4 In fact precisely this situation applies to the production of a charged Higgs boson with Mr — Mw 
e + e~ —tZ*^ W^HT , since in most models Th *C ^w- The observation of such a process would 



be a signature of an exotic Higgs sector, see for example Ref. |23 . 

5 We omit here a trivial overall branching ratio factor for the particular // final state under 
consideration. 
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Recall that as a direct consequence of the dominance of S'-wave W + W production 
at (3 < lfl 

(T ^ const, x (3 + 0{(3 3 ). (9) 
It is worth mentioning that the higher-order terms in the expansion (0) lead to a net 



negative correction to the leading ~ {3 behaviour, see for example Fig. 1 of Ref. [IS 
which originates in the high-energy SU(2)xU(l) gauge invariance cancellation. 

For E 3> T w and in the dominant ('peak') region specified by \si — m v/ \ < M W T W , 
one finds 

5 C * n - 2 arctan ( ■ (10) 

Eq. ( |I0| ) reveals the strong dependence of the coefficient 5c on the W boson virtuality 
Sx- This follows from the general nature of the Coulomb forces between unstable 
heavy particles [TT |. Thus in the large invariant mass tail (s% > Myy), 5c is strongly 



suppressed, while in the small invariant mass tail (si < Myy), 5c ~ 27r, i.e. twice 
the first-order on-mass-shell value. After integration over s\, the arctan modification 
of 5c averages to zero and the stable W result obtains. Note that while for the 
total e + e~ — > W + W~ cross section the stable W result for the Coulomb correction 
is only strongly modified by instability effects in the narrow energy region E ^ T\y 
(i.e. close to threshold), for the invariant mass distribution the arctan modification 
of 5c is essentially independent of energy for E > T w . However, far above threshold 
additional energy dependence will appear due to the screening role of the other QED 
final-state interaction mechanisms. 

For purposes of illustration, it is convenient to rewrite the cross section formula @ 
for Tw *C E <ti My/ in the dominant Si (peak) region in terms of the dimensionless 
variable x where 



x=°-± r pt~0(l). (11) 



The differential cross section in Eq. (|3p then becomes 
do ctq 1 



dx 7C 1 + x 



2 



a \ a tx 

—5h + tt arctanx 

it J p \2 



with 



E-\xT 



(12) 



(13) 



' w 

The invariant mass distribution da/dx deviates from the Breit-Wigner (1 + a; 2 ) -1 
form, corresponding to an individual W decay, because of (i) the strong dependence of 



6 Throughout this paper we assume that the ^-exchange contribution, which dominates the 
e + e~ — > W + W~ threshold cross section, is not suppressed by a particular choice of e ± beam polar- 
izations. The mass distortion effects discussed here are strictly only valid for unpolarized scattering. 
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the threshold Born cross section <t on the W momentum, and (ii) the characteristic 
behaviour of the Coulomb term discussed above. In particular, when (3 <C 1, which 
corresponds to x ~ 2E/Tw, the invariant mass distribution is strongly suppressed 
by phase space effects, Eq. (J9|). These lead to a decrease in the average value of the 
invariant mass, M = (-v/sT), in the threshold region by[] 



A Mb = at^>r w , 
AM B = 0{Y W ) at \E\ < T w , 

AM B = 0{\E\) at -E^>T W . (14) 

With increasing (3 the negative higher-order (in (3 2 ) terms in the expansion (^) become 
more and more important and, as a result, AM B changes sign (see Section |3] below). 

The characteristic dependence on the W virtuality of the Coulomb correction for 
E ^> rv always causes a decrease in the average mass compared to the Born predic- 
tion: 

AM C = 0{naT w ). (15) 

Note that this mainly arises from the Coulomb-induced asymmetry in the tails of 
the distribution, as discussed above. The shift in the actual position of the peak is 
numerically rather small. 

It is important to emphasize the difference between the predictions of Eqs. (|6l,|lO"D 
and the 'on-mass- shell Coulomb' correction, air/ {2(3). The latter would induce a shift 
in the s\ distribution towards larger values, throughout the threshold region. 

Note that there appears to be a range of collider energies around y/s ~ 190 GeV 
where the phase space and Coulomb induced distortions are of the same order in their 
effect on the average mass (see Figs. 3 and 4 below). However at this energy, which 
is of practical importance for LEP2, other mechanisms (e.g. intermediate-final and 
final-final radiative interference involving the decay products of the two W bosons) 
have to be taken into account. 

Finally, higher-order {0{a n ), n > 2) Coulomb effects could be numerically more 
important for the invariant mass distribution than for the total cross section. In 
principle, it is straightforward to take these into account using the general formalism 



presented in Ref. |TT[ (see also [f25], pb[). 

In summary, we have investigated the qualitative effects on the invariant mass 
distribution of phase-space and Coulomb corrections using a simple model in which 
only one W boson is off-mass-shell. In the following section, we shall study numerically 
the more realistic case of both W bosons being off-shell. As we shall see, the conclusions 
obtained from our model are unchanged by the more complete analysis. 



7 The subscript t B' denotes a shift in the average mass due to the Born cross section behaviour. 
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3 Invariant mass distributions in e + e — > W + W : 
quantitative discussion 

In the realistic case when both W bosons are off-shell, the W + W~ cross section can 
be written 



18} 



dsi ds 2 p(s 1 )p(s 2 ) o- (s,si,s 2 ) 

o Jo 



(16) 



where 5c is again given by Eq. and (3 = 2p/ y/s, but now with [|Tl 



V 

K 

E 



1 - 



2s(s 1 + s 2 ) - (s 1 - s 2 ) 



-M w (E + iT w ) 



AM, 



w 



(17) 



Note that we have omitted the 'hard' radiative corrections and used the 'relativistic' 
forms for p 2 and E. In what follows we will use Eqs. (|16|JT7|) to study (i) the invariant 
mass distribution da/ds\ and (ii) the average invariant mass M (which has certain 
practical advantages as an estimator of Mw |16||) defined by 



ds\ 

a(s) Jo jo 



ds 2 - (-v/si + v^s) P(si)p(s 2 ) a (s, si, s 2 ) 



1 + Tp Sc 



Figure 1 shows the normalized distribution 1/a da/ds\ (the W 'line shape') as a 
function of y/si, at three different collider energies, y/s = 165, 175 and 185 GeV.Q 
Also shown, for comparison, is the 'pure Breit-Wigner' form p(si) which corresponds, 
formally, to the y/s — > oo limit. As anticipated in the previous section, the phase space 
effects give a significant distortion to the distribution, especially close to the W + W~ 
threshold. In particular, the distribution is strongly suppressed for y/si > \/s — M w . 

The distributions in Fig. 1 include the first-order Coulomb correction 5c- To see 
the effect of this, we show in Fig. 2 the ratio 



1 da 



p(si) a ds\ 



(19) 



at sfs = 175 GeV with and without the Coulomb correction. Again we confirm 
the qualitative behaviour obtained in the model analysis of the previous section: the 



8 We use M w = 80.41 GeV/c 2 |J and T w = 2.092 GeV in our numerical calculations. All other 
parameters coincide with those used in Ref. [H (see Table 1 therein). 
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Coulomb contribution enhances (suppresses) the small (large) mass tail. Note, how- 
ever, that the effect is numerically much less significant at this energy than the dis- 
tortion due to phase space effects, which forces / to be very small for y/si > 95 GeV. 

The measurement of the W mass using the diect reconstruction method at LEP2 
involves fitting a measured invariant mass distribution, like that of Fig. 1, by a theo- 
retical distribution (in practice implemented in a Monte Carlo program) in which Mw 
is a free parameter. In this way, measured quantities like the position of the peak or 
the average invariant mass, both of which are crude measures of M^y, are corrected 
to the 'true' value. It should be clear from the above discussion that fitting the mea- 
sured distribution with a theoretical expression which does not include the final-state 
interaction effects will induce an error in the mass measurement. To quantify this, we 
focus our attention on the difference between the average W mass M flUf) defined by 
the s\ distribution and the input mass Myy, AM = M — My/. As in the previous 
section, AMb denotes the mass shift using the Born (off-mass-shell) cross section and 
AMq denotes the additional mass shift from including the 0(a) Coulomb correction. 

Figure 3 shows AM B as a function of the collider energy y/s. The behaviour can 
be understood from Fig. 1, and is exactly as anticipated in Section Near and below 
threshold, there is a strong phase space suppression for masses y/si > Mw, and so 
AMb < 0. Above threshold, the mass difference grows with increasing collider energy 
as more and more phase space for large invariant masses opens up. 

A problem with this calculation of AM is that the integrals over s\ and S2 receive 
contribution from arbitrarily small and large invariant masses (subject only to X) y/sl < 
y/s). In practice, events with very large or very small // invariant masses would not 
be classified as W decay events. In particular, lower cuts on the yfsl are required to 
eliminate rtort-W + W~ backgrounds.^ To make a more realistic calculation, therefore, 
we impose an additional cut, 



Note that this cutting procedure will to some extent complicate the calculation of the 
QED and QCD final-state radiative corrections, because of the reduction of the phase- 
space for final-state emission. The mass shifts AMb for <5 = 30 GeV and 10 GeV 
are shown as the dashed and dash-dotted curves respectively in Fig. 3. With this 
additional mass cut there is less dependence on y/s, since AMb — ► as S — > at fixed 
\/s. Note also that asymptotic values of AMb - > 0.46 (0.13) GeV are approached for 
5 = 30 (10) GeV, as ^/s oo.Q 

Figure 4 shows the additional mass shift AMc due to the Coulomb correction. 
Here we see that as long as the mass cut 5 is not particularly tight, the shift for 

9 The actual cut value will in practice depend on the particular final state, collider energy, etc. We 
choose two illustrative values for our numerical calculations. 




I = 



1,2. 



(20) 
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^/s £170 GeV is rather constant at 0(— 20) MeV. This is consistent with the model 
calculation of Section ||, which predicted a constant negative shift of order anTy/ for 
E ^> Tw- As can be derived from Eqs. ( |Tt3| , |TE| ) , AMc changes sign at lower energies 
and attains a maximum at threshold (E = 0), where the average W momentum is 
lower, (p) ~ y/My/^w- The actual maximum value depends on the cut parameter 5. 
Note once again that the mass shift decreases at fixed \fs as the invariant mass cut is 
tightened, i.e. AM C -> as 5 -> 0. 

4 Conclusions 

The success of the precision measurements of the W boson relies on an accurate 
theoretical knowledge of the details of the production and decay mechanisms. The 
favoured 'direct reconstruction' method of measuring My/ at LEP2 using the hadronic 
(qqqq) channel has an important caveat — the colour reconnection effects induced 
by the strong final-state interaction may obscure the separate identities of the W 
bosons and thus distort the mass determination []I6| . At the moment, these effects 
are not completely curable theoretically because of the lack of deep understanding of 
non-perturbative QCD dynamics. 

However, there are other effects - originating in purely QED radiative phenomena 
- which, in principle, prevent the final state being treated as two separate W decays. 
In this paper we have studied one example of this, the Coulomb interaction between 
two unstable W bosons which induces non-fact or izable corrections to the final-state 
mass distributions. Of course there is no reason why all such effects cannot, in prin- 
ciple, be computed to arbitrary accuracy in QED perturbation theory, and taken into 
account in the mass determination. In this paper we have demonstrated explicitly 
that their emission could lead to a O(20 MeV) shift in the measured mass.[^] This 
shift can only be reduced by imposing a rather tight invariant mass cut, which selects 
only those events near the peak if the distribution where the distortion is minimized. 
In particular, we have investigated the effect on the mass distribution of the QED in- 
terconnection effects generated by the first-order Coulomb corrections in the threshold 
region at LEP2. At the highest LEP2 energies, it is likely that we are overestimating 
the mass distortion effect (see for example Ref. ||). In this region, therefore, our 
results should be regarded as only a starting point for futher, more detailed studies. 
Particularly important in this respect are the QED interactions involving the decay 
products of the two W bosons, which become essential in the relativistic region. 

Finally, we note that similar non-factorizable QED final-state interaction effects 
could also be important in precision My/ measurements at the Tevatron pp collider, 

11 It is also worth mentioning that a similar effect to that described in this study could also be 
induced by final-state 'new physics' interactions, for example the exchange between the two W bosons 
of a new light scalar with a sufficiently large coupling. 
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for example in the process qg — > W(— > lu) +q. The distortion would then be manifest, 
for example, in the transverse momentum distribution of the final-state lepton. 
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Figure Captions 



[1] The distribution I /a da/dsi in e + e — > W + W — > 4/ production at y/s = 165, 
175 and 185 GeV. Also shown (dotted line) is the asymptotic form, p(si), given 
in Eq. (§. 

[2] The ratio of the mass distribution of Fig. 1 to p(si) at y/s = 175 GeV, with 
(solid curve) and without (dashed curve) the first-order Coulomb correction. 

[3] The difference between the average mass (y/si) and M\y, as a function of the 
collider energy y/s (solid curve). Also shown are the mass differences when an 
additional cut \y/sl — Mw\ < 8 (8 = 30, 10 GeV) is imposed. 

[4] The additional mass shift from including the first-order Coulomb correction, 
as a function of the collider energy y/s (solid curve). Also shown are the mass 
differences when an additional cut \y/sl — Mw\ < $ = 30, 10 GeV) is imposed. 
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